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We address the problem of noise regression in the output of gravitational-wave (GW) interferom¬ 
eters, using data from the physical environmental monitors (PEM). The objective of the regression 
analysis is to predict environmental noise in the gravitational-wave channel from the PEM mea¬ 
surements. One of the most promising regression method is based on the construction of Wiener- 
Kolmogorov filters. Using this method, the seismic noise cancellation from the LIGO GW channel 
has already been performed. In the presented approach the Wiener-Kolmogorov method has been 
extended, incorporating banks of Wiener Liters in the time-frequency domain, multi-channel anal¬ 
ysis and regulation schemes, which greatly enhance the versatility of the regression analysis. Also 
we presents the first results on regression of the bi-coherent noise in the LIGO data. 

PACS numbers: 04.80.Nn, 95.55.Ym 


I. INTRODUCTION 

Laser Interferometer Gravitational-Wave Observatory 
(LIGO) is a large-scale physics experiment targeting the 
first direct detection and study of gravitational waves 
from astrophysical sources. Two LIGO detectors in Liv¬ 
ingston, LA and Hanford, WA are being upgraded to 
increase their sensitivity by an order of magnitude and 
plan to start taking data in 2015. 

The LIGO detectors employ a Michaelson type inter¬ 
ferometer, which measures the differential length change 
induced by a passing gravitational wave in its two perpen¬ 
dicular arms with the length of 4km each pQ. While the 
fundamental sources of noise (thermal, quantum, shot, 
etc) determine the baseline LIGO sensitivity, the actual 
detector noise may also have a significant contribution 
from the environment. Numerous sources, such as ambi¬ 
ent magnetic field, power grid mains, acoustic and seismic 
disturbances are coupled into the detector sub-systems 
and can be observed as spectral lines or relatively broad¬ 
band humps in the detector output, degrading the overall 
detector sensitivity. 

In addition to the LIGO’s GW channel, data from 
thousands of Physical Environmental Monitors (PEM) 
is collected to characterize the coupling between the GW 
channel and the environment. Many environmental dis¬ 
turbances recorded in the detectors are due to the linear 
coupling to the environment, for example, the power grid 
harmonics (power lines). Also the detector data contains 
the up-conversion noise, which is produced by the in¬ 
terference of two or more environmental noises. Most 
prominent example of the up-conversion noise is the bi¬ 
linear coupling of the seismic noise which appears as side¬ 


bands around the power lines and calibration lines. To 
improve the data quality and increase the detection sen¬ 
sitivity, several methods have been used to identify and 
ultimately remove the environmental noise from the de¬ 
tector data [2'. The examples include the removal of 
seismic noise using the array of seismometers 0 , re¬ 
moval of the line noise with the finite impulse response 
(FIR) filters [4], removal of the cross-talk between an 
auxiliary channel and the channel of interest by using 
the frequency-domain linear regression [5], and predic¬ 
tion of GW strain data channel using Michelson channel 
noise [6]. 

In this paper we address the problem of the noise re¬ 
gression (prediction and cancellation) by constructing a 
special bank of the Wiener-Kolmogorov filters. The en¬ 
vironmental noise contribution to the GW channel (tar¬ 
get channel) is predicted by filtering data from one or 
more PEM channels (witness channels) simultaneously. 
To capture the frequency dependent correlation between 
the target and the witness channels and for better esti¬ 
mation of the filter parameters, the regression analysis is 
performed in the time-frequency (wavelet) domain. In¬ 
troduction of multiple witness channels in predicting the 
target channel noise increases the effectiveness of the con¬ 
structed filter banks. The strength of the presented re¬ 
gression analysis lies not only in predicting the linear, 
but also the up-conversion environmental noise. It is 
achieved by creating the synthetic witness channels from 
PEM and simulated channels, which mimic the physi¬ 
cal process of the up-conversion. This becomes useful as 
there are no witness channels that measure bilinear noise 
directly. The proposed regression model is expected to 
be more applicable when a new generation of advanced 
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detectors become operable. 

The paper is organized as follows: Section [TT] discusses 
the regression analysis; Section EH presents some appli¬ 
cations of the regression analysis; paper is concluded in 
section m 


II. REGRESSION ANALYSIS 

In application to the regression of LIGO data, we con¬ 
sider the Wiener-Kolmogorov (WK) filters [7]. Given 
sampled data from a selected auxiliary channel x mea¬ 
suring environmental noise (witness channel), a finite im¬ 
pulse response (FIR) filter a can be constructed to predict 
the noise contribution n into the GW channel t (target 
channel). The output of the filter is given by the expres¬ 
sion 


L 

n \i\ = Y a i x i+h (2- 1 ) 

j=~L 

where aj are the filter coefficients and 2L + 1 is the filter 
length. The filter is obtained by minimizing the mean- 
square error 


N+L / L 

x 2 = Y 'W “ J2 

i=L \ j=—L 

where N is the number of data samples used for the es¬ 
timation (training) of the filter. The filter coefficients aj 
are calculated by solving the Wiener-Hopf (WH) equa¬ 
tions 



Rxx& — P tx- (2*3) 

where p tx is the vector with 2L + 1 components repre¬ 
senting the cross-correlations between the witness chan¬ 
nel with the target channel and R xx is (2L + 1) x (2L + 1) 
auto-correlation matrix for the witness channel. The ma¬ 
trix R xx is positive definite and therefore non-singular, 
yielding a unique solution for the filter coefficients. Usu¬ 
ally, the solution is obtained with the Levinson-Durbin 
algorithm [8] without the explicit inversion of R xx . The 
predicted noise is calculated by using equation |2.1| and 
can be subtracted from the target channel, thus reducing 
the noise in the GW data. 

A special case of the Wiener filter is the linear pre¬ 
diction error filter, which is obtained by minimizing the 
following equation 

N+L / L \ 2 

x 2 = Y C N+j ■ (2.4) 

i=L y j=-L,j^ 0 / 

In this case, the data sample i = L is predicted by the 
surrounding samples in the same time series. This is well 
suited for the prediction of a quasi-monochromatic noise 
in the target channel. 


The RMS value of the witness channel data can be 
broken down into coupled and uncoupled parts. The tar¬ 
get and witness channel data are whitened before being 
presented to the analysis, e.g., has a RMS value of 1. 
Following this, one can write, 

i = r l+ r li ( 2 - 5 ) 

where r c /r u is the RMS of the coupled/uncoupled part 
in the witness data. 


A. Regression with WDM 

The interferometer data span a wide frequency band 
(0-8 kHz) and have a large dynamic range. Therefore, 
construction of the WK filters in the time domain re¬ 
quires long filters. For example, to regress power line har¬ 
monics (15 in the 0-1 kHz band) a WK filter with approx¬ 
imately 10000 coefficients should be constructed. Apart 
from the computational complexity associated with the 
inversion of the 10000x10000 matrix, one need to de¬ 
termine accurately 10000 filter coefficients, while only 
30 parameters (amplitudes and phases of 15 harmon¬ 
ics) are relevant. Also, the WK filter is affected by the 
spectral leakage, which may fail to capture all details of 
the correlation between the witness and the target chan¬ 
nels. To solve these problems, we propose to use the fast 
Wilson-Daubechies-Meyer transformation (WDM) [9] to 
split data into the frequency sub-bands (time-frequency 
map). Figure [l] shows an example of the WDM time- 
frequency map for a segment of data from the Hanford 
detector. The WDM transformation is orthonormal, in¬ 
vertible and has a low spectral leakage, which makes it 
a unique tool for the regression analysis. The data in 
each frequency sub-band (WDM layer) are a time se¬ 
ries used for the construction of a single WK filter as 
described in the beginning of this section. Therefore, in¬ 
stead of one long WK filter, a bank of much shorter WK 
filters is constructed. Each filter is trained individually 
to capture details of the target-witness correlation in the 
corresponding sub-band with just a few filter coefficients. 
For example, a typical analysis in the frequency band 0- 
2048 Hz involves the construction of 2048 independent 
WK filters with 11 coefficients each. Therefore, instead 
of the inversion of a very large matrix, one has to invert 
2048 much smaller matrices, which significantly reduces 
the complexity of the regression problem. By applying 
the inverse WDM transform, clean data h — x in the time 
domain can be obtained. 


B. Multiple witness channels 

The GW channel gets contaminated by the environ¬ 
mental noise entering from different physical locations 
and directions. It is a commonly encountered scenario 
when a particular environmental noise coupled to the tar¬ 
get channel is measured by the witness channels situated 
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FIG. 1: An example of time-frequency map obtained by ap¬ 
plying WDM transform on LIGO data (not all layers have 
been shown for clarity purposes). Layers have a bandwidth 
of 1 Hz. Regression analysis is performed for each layer sepa¬ 
rately. Once the prediction is calculated and subtracted from 
each frequency sub-band of the target data, the inverse trans¬ 
form is applied to bring the data back to the time domain. 


at different locations. In general, when making a predic¬ 
tion over a period of time, data from multiple witness 
channels are expected to provide more complete infor¬ 
mation about the environment and improve the estima¬ 
tion of the filter coefficients. For example, the seismome¬ 
ters are installed at various key locations in the groups 
of three throughout the LIGO sites, measuring seismic 
noise along specified x, y and z directions. Therefore, a 
prediction of the seismic contribution to the target chan¬ 
nel would require a simultaneous use of all three seismic 
channels. The regression analysis addresses the multiple 
witness channel case by extending Equation |2.2| 

N+L / L L 

x 2 = E (- E a i Xi +i - E bkVi + k 

i=L \ j= — L k= — L 

L 

^ ^ CmZi-\-m T 

m=—L 



Where h is the target channel, x, y, z are the witness 
channels and a, b, c are the filter coefficients. The corre¬ 
sponding Wiener-Hopf equation is 
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(2.7) 


where R xx /C xy are Hermitian Toeplitz matrices of the 
witness auto/cross correlation and p^, p hy , Ph z are 


Number of channels 

1 

2 

4 

8 

16 

RMS of prediction 

.07 

0.09 

0.13 

0.19 

0.26 


TABLE I: In the absence of correlation between the target and 
the witness channels the RMS of the prediction increases with 
the number of witness channels due to over-fitting, i.e. the 
analysis makes prediction even in the absence of correlation. 


cross-correlation vectors between the witness and the tar¬ 
get channels. The output of the multi-channel filter is 
given by: 

L L L 

x\i\ = y ^ CLjXi+j H- y ^ bhi/i+k T y ^ GnTi+m- 

j=—L k=—L m=—L 

( 2 . 8 ) 

The example of improvement in prediction due to the 
multi-channel regression is presented in figure |2j It shows 
the regression of seismic noise using the “coil current” 
(CC) channels. The CC channels measure the current in 
the electromagnets driving small magnets attached to the 
test masses. The main function of these electromagnets 
is to counteract the seismic motion and hold the test 
masses in place. Since seismic noise at the sites usually 
affects nearby test masses, at a given time, only a small 
subset out of total 16 CC channels (positioned at the 
initial and end test mass) exhibits correlation with the 
GW channel. Therefore, the regression with a single CC 
channel is not effective. One can construct 16 regression 
cases and apply them consecutively, one CC channel at a 
time. However, as figure [2] shows, the best performance 
is obtained with the multi-channel regression when all 16 
channels are used to construct a single regression filter. 

The multi-channel analysis |2.6| presents two shortcom¬ 
ings. First, if the witness channels are highly correlated, 
the matrix Q can be rank deficient. Second, at any given 
moment a significant fraction of witness channels may not 
have any measurable correlation with the target channel 
and just add noise at the filter output. This is demon¬ 
strated in table |TJ where we consider 16 witness channels 
with random Gaussian noise. In this case the output of 
the filter is also Gaussian noise with the RMS increas¬ 
ing as y/~M, where M is the number of the noisy witness 
channels used in the regression. This case shows, that for 
a significant number of noisy witness channels the regres¬ 
sion analysis may result in the over-fitting of the target 
data. 


C. Regulators 

Both shortcomings can be addressed by using regula¬ 
tors. The equation |2.7| can be rewritten as 
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FIG. 2: Prediction of seismic noise using 16 coil current chan¬ 
nel. Black curve is the original power spectrum. Red curve 
is obtained by conditioning the data using multichannel anal¬ 
ysis. Blue curve is obtained by performing single channel 
regression analysis, on the target data, using the CC chan¬ 
nels consecutively. Purple curve is a single channel analysis 
with a different order of the CC channels. 


where O is the orthogonal matrix that diagonalizes Q, 
and Xi (i = 0,1,2...) are the eigenvalues arranged in 
the decreasing order. The eigenvalue distribution cap¬ 
tures the principal regression components and identifies 
a strong correlation in the witness data. Figure [3] plots 
the eigenvalue distributions for different regression cases. 
Because the witness data are whitened before performing 



FIG. 3: Eigenvalue distribution of matrix Q for some cases. 
Black - all 16 witness channels carry the same harmonic noise 
along with additional white noise. Red - all the 16 witness 
channels carry white noise, none of the witness channels are 
correlated with each other. Green - all the 16 witness channels 
carry white noise and all are correlated with each other. 

the analysis, the mean of the eigenvalue distribution is 


always equal to one. When A i » 1, this is an indication 
of strong narrow-band noise artifacts in the witness data. 
Respectively, when A i « 1, the witness data are highly 
correlated and the matrix Q has a lower rank. The small 
eigenvalues should be regulated to avoid un-physical so¬ 
lutions of the WH equation. The regulated solution is 
obtained by modifying the matrix A, 
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( 2 . 10 ) 

where the eigenvalues below the threshold X t a set 
to some constant C. We distinguish three regulation 
schemes: hard ( 1/C = 0), mild (C = Ao) and soft 
(C = A p ). The proposed regulators handle the problem 
of the rank deficiency in the matrix Q by suppressing 
the eigenvectors with small eigenvalues thereby prevent¬ 
ing un-physical solutions. In the case of no correlation 
between the witness and the target channels, the regula¬ 
tors also constrain the system to reduce the RMS of the 
random noise in the output of the filter. 


D. Effect of regulators on prediction 


The regression parameters and the type of the regula¬ 
tor to use depend on the type of problem on hand. The 
coupling of the environmental noise to the target and 
the witness channels usually changes with time. Hence, 
no general prescription can be made on the application 
of the regulation schemes. Below, we discuss figures of 
merit for some frequently encountered cases. In all dis¬ 
cussed cases, the target and the witness channels differ 
by an additional Gaussian noise in the witness channel. 

1) All three regulators are very close in performance 
for the prediction of quasi-monochromatic noise. Fig¬ 
ure [4] compares the regression of simulated monochro¬ 
matic noise with different regulators. As only few eigen¬ 
values are significant in this case, all regulation schemes 
remove the noise peak. The effect of the over-tuning is 
clearly visible around the peak when no regulator is used. 
The hard regulator is the best choice for such cases. 

2) In the case when all witness channels measure inde¬ 
pendent broadband noise, the eigenvalue distribution is 
Gaussian. If none of the witness channels is correlated 
with the target channel, the goal is to reduce the RMS 
value of the prediction. Table [IT] shows the performance 
of various regulators. In this case using a hard regulator 
with a high threshold works the best. On the contrary, 
if the correlation exists only between one witness chan¬ 
nel and the target channel at a time, regulating a large 
number of eigenvectors results in reduced performance, 


as shown in table III A soft regulator or no regulator 
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FIG. 4: Prediction of simulated monochromatic noise using 
different regulators with \ t = 1. Black curve is the original 
power spectrum. Red curve is obtained after conditioning 
the data. The 99.5 Hz line has been cleaned using the hard 
regulator, the 101.5 Hz line has been cleaned using the mild 
regulator, the 103.5 Hz line has been cleaned using the soft 
regulator, and the 105.5 Hz line has been cleaned with no 
regulators. 


Threshold 

0.25 

0.5 

0.75 

1.0 

1.25 

1.5 

Hard 

0.26 

0.26 

0.23 

0.18 

0.12 

0.05 

Mild 

0.26 

0.26 

0.24 

0.20 

0.18 

0.17 

Soft 

0.26 

0.26 

0.24 

0.24 

0.18 

0.18 


TABLE II: Effect of regulators when all the 16 witness chan¬ 
nel carry Gaussian noise and none of them are correlated with 
the target channel. Table shows the RMS value of the predic¬ 
tion. 

works best for this case. This case is only relevant when 
the correlation between the witness and the target chan¬ 
nel is strong. If the correlation is weak, due to the con¬ 
tamination by uncorrelated witness channels, regression 
analysis does not perform well; for instance, with 16 wit¬ 
ness channels and one of them correlated with the target 
channel with r c = .6, the RMS of the residual is 0.8. 

3) In the case when all the witness channels have Gaus¬ 
sian noise present and more than one are correlated with 
the target channel, there is splitting in the eigenvalue dis¬ 
tribution (figure [3| . All the regulators work equally well 
in this case. For instance, when 8 out of 16 witness chan¬ 
nels are correlated with the target channel with r c = .6, 
the RMS of residual is close to .4 for all the regulators 
and threshold values of up to At = 1.5. 

III. APPLICATIONS OF WDM REGRESSION 

There is no one single prescription for the WDM re¬ 
gression. Each regression case depends on a set of avail¬ 
able witness channels and type of the noise to be pre¬ 
dicted. Primary consideration is that witness channels 


Threshold 

0.25 

0.5 

0.75 

1.0 

1.25 

1.5 

Hard 

0.152 

0.152 

0.381 

0.617 

0.811 

0.953 

Mild 

0.152 

0.152 

0.262 

0.352 

0.396 

0.410 

Soft 

0.152 

0.152 

0.161 

0.210 

0.283 

0.362 


TABLE III: Effect of regulators when only one witness 
channel out of 16 is correlated with the target channel with 
r c = 0.99. All channels contain Gaussian noise. Table shows 
RMS of the residual. 



FIG. 5: Performance of various regulators (with threshold 
0.5) in subtraction of seismic noise using 16 CC channels. 
Black curve is the original power spectrum. Colored curves 
are obtained after conditioning the data. Red - No regulator, 
Green- Soft regulator, Blue- Mild regulator, and Cyan - Hard 
regulator. Effect of regulator is clearly visible. It can be 
deduced that only one CC channel is correlated with the GW 
channel at a time. 

should not be sensitive GW. This consideration can be 
lifted, when cleaning upconverted noise, if the frequency 
domain of synthetic channels (discussed later) does not 
overlap with the frequency domain of the channels used 
in their construction. Below, we consider several typical 
regression cases performed on the data obtained from the 
Hanford detector. 


A. Regression of noise with linear coupling. 

One of the most visible noise artifacts in LIGO detec¬ 
tors are power harmonics showing up in the noise spec¬ 
trum as the quasi-monochromatic lines at multiples of 
60 Hz. They can be monitored with the voltage moni¬ 
tors and magnetometers. One or several voltage monitors 
can be used as witness channels. Figure [6] shows the noise 
spectrum before and after the application of the WDM 
regression analysis. Note, that the line adjacent to 120 
Hz in figure [6] (middle) is not removed because it is not 
monitored by the voltage channels. 

Other frequently visible features in the LIGO noise 
spectrum are due to the mechanical resonances of the 
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FIG. 6: Cleaning of power-lines using voltage monitor 

(H0:PEM-LVEA2_V1) as the witness channels. Black curve 
is the original spectrum (H1:LDAS-STRAIN); red is obtained 
after subtracting the prediction. 


seismic stacks. These lines are much wider than the 
power lines, but they can still be efficiently removed by 
using regulated regression with multiple witness channels 
(Figure [7]) 



FIG. 7: Subtraction of mechanical resonances using 14 ac¬ 
celerometers and 2 microphones. Black curve is the original 
spectrum (H1:LDAS-STRAIN); red curve is obtained after 
subtracting prediction using a mild regulator (with threshold 
1 ). 


ITMX_COIL_UR, see Section IIB). We do not observe 
a significant rank deficiency of the witness matrix indi¬ 
cating that there is no significant correlation between the 
CC channels. As expected, in this case, the best regres¬ 
sion is obtained when no regulator is used. 


B. Regression of non-linear noise 

Some noise artifacts in LIGO are produced by the in¬ 
terference of two (bi-linear noise) or more noise sources 
in the detector. The majority of the bi-linear noise cases 
are due to the up-conversion of the low frequency seis¬ 
mic noise, which is observed as the side-bands around 
the power lines, calibration lines, violin modes and other 
narrow-band features in the noise spectrum. There are 
no direct witness channels in LIGO to monitor such bi¬ 
linear noise. As the CC channels provide a good measure 
of the low frequency seismic noise, they can be used along 
with the other channels to synthesize the bi-linear witness 
channels. Such synthetic witness channels are obtained 
by multiplying the CC time series by the carrier time 
series representing a narrow-band spectral artifact. For 
example, when the carrier is one of the power lines, the 
bi-linear witness channel can be obtained as a product of 
the time series from a CC channel and a voltage moni¬ 
tor or magnetometer. Figure [8] shows the example of the 
up-conversion noise regression around 180 Hz power line 
with 16 synthetic witness channels constructed from the 
CC channels and one voltage monitor channel. Usually, 
such cases need to be regulated to avoid the regression 
artifacts. 



A more difficult linear regression case is the predic¬ 
tion of the broadband noise, like seismic noise. Figure [2] 
shows the example of the seismic noise regression using 
16 coil CC channels (Hl:SUS-ITMX_COIL_LL, H1:SUS- 


ITMX_COIL_LR, 

ITMX_COIL_UR, 

ETMX_COIL_LR, 

ETMX_COIL_UR, 

ITMX_COIL_LR, 

ITMX_COIL_UR, 

ITMX_COIL_LR, 


Hl:SUS-ITMX_COIL_UL, 

Hl:SUS-ETMX_COIL_LL, 

Hl:SUS-ETMX_COIL_UL, 

Hl:SUS-ITMX_COIL_LL, 

Hl:SUS-ITMX_COIL_UL, 

Hl:SUS-ITMX_COIL_LL, 

Hl:SUS-ITMX_COIL_UL, 


H1:SUS- 

H1:SUS- 

H1:SUS- 

H1:SUS- 

H1:SUS- 

H1:SUS- 

H1:SUS- 


FIG. 8: Removal of the up-conversion noise around 180 Hz 
power line with 16 synthetic witness channels constructed 
from 16 coil current channels and the voltage monitor[6] Black 
curve is the original power spectrum. Red curve is obtained 
after conditioning the data. 

Often, the carrier time series are not readily mea¬ 
sured by any PEM channels. If the carrier frequency 
is quite stable, the carrier witness channel can be simu¬ 
lated with the monochromatic signal. Figure [9] shows how 





















7 


this method works for the regression of the up-conversion 
noise around the calibration lines. Also, the carrier wit¬ 
ness channel can be extracted from the strain (target) 
data with the linear prediction error (LPE) filter. In this 
case, the LPE filter is trained on the WDM data in the 
narrow-band 1 Hz) around the carrier frequency. The 
output of the LPE filter is used as the carrier witness 
channel. Figure [To| shows the subtraction of the seismic 
up-conversion noise around the beam jitter peaks, which 
were not monitored by any witness channels during the 
LIGO science runs. 



FIG. 9: Removal of up-conversion around 113 Hz calibration 
line using coil current channels and a simulated noise. Black 
curve is the original power spectrum. Red curve is obtained 
after conditioning the data. 



in the monitoring of the noise artifacts present in the 
target data. However, many channels do not exhibit any 
significant correlation with the strain channel or this cor¬ 
relation varies significantly over the time. The proposed 
regression analysis can be used to monitor the correla¬ 
tion between the target and the witness channels during 
the data taking runs. The measure of the correlation is 
the RMS value of the whitened prediction signal, which 
should be between 0 (no correlation) and 1 (strong corre¬ 
lation). Since it is hard to visualize all RMS values from 
hundreds of LIGO channels, the maximum RMS value 
among all the WDM layers in each witness channel can be 
used to identify the strongest correlation with the strain 
channel. Figure [Tl] shows the maximum RMS value for 
approximately four hundred LHO channels (all the chan¬ 
nels which have sampling rate of 2048 Hz or higher). To 
obtain this figure of merit, for each witness channel the 
2048 regression filters were constructed covering a 1 Hz 
frequency sub-band each. The RMS values of the filter 
output were calculated and the maximum value is dis¬ 
played as a function of the channel identification number 
and the GPS time of the test data segment. In this plot, 
the contribution of the power lines has been excluded, as 
they are correlated with many channels. There is no sig¬ 
nificant correlation for most of the tested channels: green 
and blue areas in the plot where RMS<0.5. Some aux¬ 
iliary channels have a strong correlation with the strain 
channel: yellow and red areas in the plot. 



FIG. 10: Removal of up-conversion around jitter peaks using 
coil current channels and self-predicted carrier. Black curve 
is the original power spectrum. Red curve is obtained after 
conditioning the data. 


C. Monitoring of Gravitational Wave Data 

LIGO has hundreds of auxiliary channels. Some of 
these channels (like voltage monitors) are very effective 


FIG. 11: The maximum RMS value as a function of the chan¬ 
nel identification number and GPS time. 

There are other possible figures of merit that can be 
used. For example, the average RMS value of the pre¬ 
dicted noise from several WDM layers is a good indica¬ 
tor of the broad-band correlation with the strain data. 
Also, the witness channels can be monitored in groups 
by using the multi-channel regression. This monitoring 
method may help to identify channels, which are weakly 
correlated with the GW channel individually, but may 











































reveal a strong correlation as a group. 

IV. CONCLUSION 

We present a novel method for regression of LIGO 
data. Using the LIGO PEM data estimates the contri¬ 
bution of the environmental noise into the GW chan¬ 
nel. The regression analysis is performed in the time- 
frequency domain obtained with the WDM transform. 
Banks of the Wiener-Kolmogorov filters are constructed 
to capture the details of the witness-target correlation. 
The regression analysis is extended to multiple channel 
case as, in general, a single witness channel may not have 
sufficient information about the correlation between the 
GW channel. Inclusion of multiple channels is shown to 
increase the efficiency of the regression analysis, and also 
present the problems of over-fitting and rank deficiency 
of linear equations. Regulators are introduced to miti¬ 
gate these problem, thereby, removing the need for the 
micromanagement of multiple channels used in the anal¬ 
ysis. Linear noise is regressed by directly accepting the 
data from the PEM channels as the witness channels, 
while bi-linear noise is handled by constructing witness 


channels from two or more PEM and/or simulated chan¬ 
nels. Regression in the WDM domain along with the 
freedom to adjust the regulators gives the analysis the 
flexibility in handling variety of noises appearing at vari¬ 
ous frequencies. Regression analysis can also be used for 
monitoring data from future LIGO detectors. With hun¬ 
dreds of PEM expected to be functional for future detec¬ 
tors, channel monitoring will identify the PEMs having 
the greatest effects on a detector’s performance. 
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